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ABSTRACT 

We present the elemental abundance and H2 content measurements of a Damped Lyman-a 
(DLA) system with an extremely large Hi column density, logA^(Hi) (cm~^) = 22.0±0.10, at 
Zabs = 3.287 towards the QSO SDSS J 081634+144612. We measure column densities of H2, Ci, 
C I* , Zn n, Fe 11, Cr 11, Ni 11 and Si 11 from a high signal-to- noise and high spectral resolution VLT- 
UVES spectrum. The overall mctallicity of the system is [Zn/H] — —1.10 ±0.10 relative to solar. 
Two molecular hydrogen absorption components arc seen a.t z ~ 3.28667 and 3.28742 (a velocity 
separation of ~ 52 km s~^) in rotational levels up to J = 3. We derive a total H2 column density 
of logiV(H2) (cm-2) = 18.66 and a mean molecular fraction of / 2iV(H2)/[27V(H2) + iV(Hi)] = 
^0-3.04±o.37^ typical of known H2-b caring DLA systems. From the observed abundance ratios we 
conclude that dust is present in the Interstellar Medium (ISM) of this galaxy, with a enhanced 
abundance in the H2-bearing clouds. However, the total amount of dust along the line of sight is 
not large and does not produce any significant reddening of the background QSO. The physical 
conditions in the H2-bearing clouds are constrained directly from the column densities of H2 in 
different rotational levels, C i and C i"^. The kinetic temperatme is found to be T « 75 K and the 
particle density lies in the range tih = 50—80 cm^'^. The neutral hydrogen column density of this 
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DLA is similar to the mean H i column density of DLAs observed at the redshift of 7-ray bursts 
(GRBs). We explore the relationship between GRB-DLAs and high column density end of QSO- 
DLAs finding that the properties (metallicity and depletion) of DLAs with log (Hi) > 21.5 in 
the two populations do not appear to be significantly different. 

Subject headings: quasars: general — quasars: absorption lines — ISM: molecules 



1. Introduction 

Despite accounting for only a small fraction of 
all the baryons in the Universe (see, e.g., Petit- 
jean et al. 1993), the physical state of the neutral 
and molecular phases of the interstellar medium 
is a crucial ingredient of galaxy formation. These 
gaseous phases are at any redshift the reservoir of 
gas available for star formation. At high redshift, 
most of the neutral hydrogen mass is revealed by 
the damped Lyman-a absorption system (DLAs) 
detected in the spectra of background quasars (see 
e.g. Wolfe et al. 2005 for a review). Since DLAs 
are easy to identify in QSO spectra and the Hi 
column densities can be measured accurately, it 
is possible to derive the cosmological mass den- 
sity of the neutral gas at different redshifts, in- 
dependent of the exact nature of the absorbers, 
provided a sufficiently large number of background 
quasars are observed (see Prochaska & Wolfe 2009, 
Guimaraes et al. 2009, Noterdaeme et al. 2009b). 

Key results of DLA surveys include the indi- 
cation that the Af(Hi) frequency distribution de- 
viates significantly from a single power-law with 
/(A(H i)) sharply steepening at log iV(H i) > 21.5. 
However, Zwaan & Prochaska (2006) used CO 
emission maps in the nearby Universe to show 
that the H2 column density distribution function 
is a continuous extension of the H i distribution 
for high column densities. The transition happens 
at log A(Hi) ^ 22 which is the approximate col- 
umn density associated with the conversion from 
Hi into H2 (e.g., Schaye 2001). 

The slope of the column density distribution 
at large 7V(Hi), a ^ 3.5, implies that systems 
with very large column density are very rare. In- 
deed, until very recentl}0, only one DLA system 
with logA^(Hi) ~ 22 was reported in the liter- 
ature: the system at Zabs = 3.287 towards SDSS 
J 081634-f 144612 (from the SDSS-II DLA catalog. 



^In the last stages of this work, the discovery of another 
DLA with log A'^(H l) > 22 was reported by Noterdaeme et 
al. (2012) and Kulkarni et al. (2012). 



Noterdaeme et al. 2009b). However, DLAs with 
such high column densities are frequently detected 
at the redshift of 7-ray bursts (e.g. Savaglio et al. 
2003; Vreeswijk et al. 2004; Fynbo et al. 2006; 
Jakobsson et al. 2006; Prochaska et al. 2007; 
Fynbo et al. 2009; Ledoux et al. 2009; Savaglio 
2010). Therefore, comparing the chemical and 
physical properties of the strongest QSO-DLAs to 
that of GRB-DLAs may provide clues to under- 
stand the nature of these absorbers. 

We present here a detailed study of the gas- 
phase abundances of metals, dust and molecules 
in the QSO-DLA towards SDSS J 081634-f 144612, 
based on high spectral resolution data. The pa- 
per is organized as follows. In Sect. 2 we pro- 
vide details of observations and data reduction. 
The hydrogen column density of the DLA, the 
metal, dust and molecular content are discussed in 
Sect. 3. The physical state of the H2 bearing gas 
component is discussed in Sect. 4. In Sect. 5, we 
finally investigate the connection between QSO- 
DLAs and GRB-DLAs. 

2. Observations 

We observed quasar SDSS J 081634-f 144612 
twice, in September 2008 and April 2009, with 
the high-resolution Ultraviolet and Visual Echelle 
spectrograph (UVES, Ballester et al. 2000) 
mounted on the ESO Kueyen VLT-UT2 8.2 m tele- 
scope at Cerro Paranal, Chile. Observations have 
been performed under programs 081.A-0334(A), 
PI. S. Lopez in visitor mode, and 282.A-5030(A), 
PI. P. Noterdaeme in service mode. Ten expo- 
sures were taken for a total of 12.4 hours exposure 
time: nine exposures using Dichroic 2 with a set- 
ting 437-1-760 nm plus one 5400 s exposure with 
the red arm centered at 550 nm that covers the 
Lyman-a absorption. A slit width of 1 arcsec and 
2x2 pixel binning were used, resulting in a spectral 
resolution of 50,000. 

The quasar spectrum was reduced using the 
UVES pipeline (see e.g. Ledoux et al. 2003 for 
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details). The main characteristics of the pipeUne 
are to perform a precise inter-order background 
subtraction, especially for master flat-fields, and 
to allow for an optimal extraction of the object 
signal rejecting cosmic rays and performing sky 
subtraction at the same time. The pipeline prod- 
ucts were checked step by step. The wavelength 
scale of each reduced spectrum was then converted 
to vacuum-heliocentric values and the spectra re- 
binned to a constant wavelength step. No further 
rebinning was performed during the analysis of the 
whole spectrum. Individual 1-D exposures were 
scaled, weighted and combined together. 

In order to derive the physical parameters of 
the absorption features, we fit the metal absorp- 
tion profiles with multiple Voigt profiles using 
VPFIT (Carswell et al. 1987). The continuum 
level was obtained locally in the vicinity of each 
metal absorption feature. Molecular hydrogen fea- 
tures were fitted altogether using fit/lyman and 
after normalizing the corresponding region of the 
Lyman-a forest. Atomic data for metal species 
and H2 are, respectively, from Morton (2003) and 
Bailly et al. (2010). In the following, solar abun- 
dances are taken from Lodders (2003). The origin 
of the velocity scale {v = km s~^) is set at the 
redshift of the single Ci component, z = 3.28746. 

3. Abundances 

3.1. Hi and metal content 

The neutral hydrogen column density was mea- 
sured from the fit of the damping wings of the 
Lyman-a absorption at z ~ 3.287. We find 
logA^(Hi) (cm-2) = 22.0±0.10. The observed 
Damped Lyman-a absorption together with the 
best fitted Voigt profile is shown in Figure [TJ The 
dashed lines indicate the profiles corresponding to 
logA^(Hi) = 21.9 and 22.1. 

We detect absorption lines of Ci, Ci*, Znii, 
Pen, Cm, Niii and Si 11, spread over about 
150 kms~^. As can be seen in Fig. [2l the ab- 
sorption profiles are not strongly saturated except 
may be for SiiiAlSOS. We are therefore confi- 
dent that our column density determinations are 
robust. The fit to the absorption lines are over- 
plotted on the Figure as red solid lines. The singly 
ionized species are expected to be the dominant 
contributors to the abundances of the correspond- 
ing elements in Hi clouds with such high column 



densities. 

Voigt profile fitting is performed simultane- 
ously for all the absorption lines keeping the same 
number of components having same redshifts and 
Dopplcr parameters for all singly ionized species. 
Note that ZniiA2026 is blended with MgiA2026. 
The contributions of the latter is taken into ac- 
count in the fits and found to have negligible in- 
fluence on the derived N(Znii). 

Eight velocity components are necessary to 
model the proflles (see Fig. [2]). We report in Ta- 
ble [1] the results of the fits, column density and 
Doppler parameter, for each of the components. 
Because of saturation effects, the Sin column den- 
sity should be considered a lower limit. However, 
given the shape of the absorption, we are confi- 
dent that the true value cannot be much larger. 
The column density of singly ionized iron, despite 
their blending, could be reliably measured. 

Although the profile decomposition may not 
be unique, three distinct clumps can be identified 
(#1, 2 and 3 from blue to red) at mean redshifts of 
3.28661, 3.28735 and 3.28814, and made of respec- 
tively, 3, 2 and 3 components. H2 absorption is de- 
tected in ^1, which is the strongest clump in metal 
species, but most of the H2 is found associated to 
the C I component in clump ^2. Unfortunately, 
we cannot determine iV(Hi) in each clump. The 
mean metallicit}!! in the cloud, derived by adding 
the column densities of the eight components are: 
[Zn/H] = -1.10 ± 0.10, [Si/H] > -1.23 ± 0.10, 
[Cr/H] = -1.58±0.10, [Ni/H] = -1.68±0.10 and 
[Fe/H] = -1.58 ±0.10. 

3.2. Dust 

In the ISM of the Galaxy, zinc is virtually unde- 
pleted onto dust grains when Si, Cr and Ni are. We 
find in the present DLA mean relative abundances: 
[Si/Zn] = -0.13 ± 0.03, [Cr/Zn] = -0.48 ± 0.02, 
[Ni/Zn] = -0.58±0.02and [Fe/Zn] = -0.48±0.02 
indicating that the overall depletion of Si, Cr, Ni 
and Fe is similar to what is seen in the gas from 
the halo of our Galaxy (Welty et al. 1999). This is 
also typical of DLAs where H2 is detected (Noter- 
daeme et al. 2008). From Fig. [31 which presents 
the depletion patterns relative to zinc observed 
component by component, it is clear however that 



■^Metallicities are given relative to solar: [X/H] = 
log(7Vx/A'm)DLA - log{X/H)Q 
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the depletion is enhanced in the main H2-bearing 
component. This situation is similar to that of the 
H2-bcaring component towards Q 0013—004 (Pe- 
titjean et al. 2002), where higher depletion factors 
are seen in the H2 components. 

From the flux-calibratcd SDSS spectrum, it is 
also possible to estimate the reddening induced 
by the presence of dust in the DLA to the back- 
ground QSO light, following the method described 
in Noterdaeme et al. (2009b). In Fig. H we 
show that the SDSS spectrum of J 081634+144612 
is well matched with the SDSS composite spec- 
trum from Vandcn Berk et al. (2001), shifted 
to the same emission redshift and reddened us- 
ing a SMC extinction law (Gordon et al. 2003) 
at z = 3.286 with E(B-V) = 0.05±0.06. The as- 
sociated uncertainty is obtained from the disper- 
sion measured for a control sample of 163 SDSS 
QSOs from Schneider et al. (2010) with emission 
redshift within ±0.02 to that of J 081634-^144612. 
This means that there is no significant redden- 
ing of the quasar J 081634+144612. Overall, the 
measured extinction-to-dust ratio, Av/A^(Hi) < 
5 X 10^^"^ magcm^ (2 cr) is typical of that of the 
general DLA population (Vladilo et al. 2008). 
This, together with the presence of H2 in a com- 
ponent with higher depiction factor suggests that 
appreciable fraction of H i may be associated with 
components that do not have H2. 

3.3. H2 content 

Molecular hydrogen is detected in two distinct 
sub-systems at Zabs = 3.28667 and 3.28742, sepa- 
rated by ~52 km with absorption lines from 
rotational levels up to J = 3 (see Figure [S]). The 
results of the fits to the numerous absorption lines 
are given in Table [2] The Zabs = 3.28742 compo- 
nent alone contains about 90% of the total H2 col- 
umn density in the absorber, log iV(H2) = 18.62, 
and coincides with the Ci component. This is 
expected because the energy of the photons that 
ionize C i is close to that of photons that dissoci- 
ate H2. The total H2 column density integrated 
over the two components and all rotational lev- 
els is log iV(H2) = 18.66 ± 0.27 (cm'^) corre- 
sponding to a molecular fraction of log / = log 2 x 
iV(H2)/(2x iV(H2) + A^(Hi)) = -3.04 + 0.37 if we 
assume that the totality of neutral hydrogen is as- 
sociated with the H2 components. This value is 
amongst the lowest observed in H2 bearing DLAs 



with metallicities [Zn/H] > —1.3 (Pctitjean et al. 
2006, Noterdaeme et al. 2008). This again may 
indicate that appreciable fraction of H I may be 
associated with components that do not have H2 
(see also Noterdaeme et al. 2010, Srianand et al. 
2010, 2012). 

4. Physical state of the gas 

4.1. Excitation of H2 

From the detection of H2 in different rotational 
levels (J=0 to J=3, see Table it is possible 
to put constraints on the physical state of the 
gas. The excitation temperature Tqj between ro- 
tational levels and J is defined as 

N{J) _g{J) 

-B(OJ)/kToj /-|^ 

N{0) 5(0) ^ ^ 

where g{3) is the statistical weight of the rota- 
tional level J: ^(J) = (2J + 1)(2I+ 1) with nuclear 
spin I = for even J (para-H2) and I = 1 for 
odd J (ortho-H2), k is the Boltzmann constant, 
and E{OJ) is the energy difference between level J 
and the ground state (J=0). If the excitation pro- 
cesses are dominated by collisions, then the popu- 
lations of the rotational levels follow a Boltzmann 
distribution described by a unique excitation tem- 
perature for all rotational levels. This is gener- 
ally the case for low rotational levels which have 
a de-excitation time scale larger than the collision 
time-scale. Indeed, Tqi is a good indicator of the 
kinetic temperature (Roy et al. 2006; Le Petit et 
al. 2006) especially in clouds similar to the one we 
study here (log A^(Il2) > 18). However, because of 
the small energy difference between the J=:0 and 
J=l levels, the value of Tqi is very sensitive to 
uncertainties on A^(Il2,J=0) and A^(H2,J=1) and 
the use of higher rotational levels may help derive 
a better constraint on Tk. From Fig. [6j it can be 
seen that the population of J=0 to J=2 levels can 
be described by a unique excitation temperature 
Tex = 69tl° K and Tex = 79tlt K for the first and 
second component, respectively. 

These temperatures are slightly smaller than 
what was found in previous studies of Il2-bearing 
DLAs (e.g. Ledoux et al. 2003, T - 90 to 180 K; 
Srianand et al. 2005, T - 153+78 K), but simi- 
lar to what is measured in the ISM of our Galaxy 
(77+17 K; Rachford et al. 2002) and in the Mag- 
ehanic Clouds (82+21 K; Tumhnson et al. 2002), 
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where Hi column densities are also large. Note 
that, temperatures observed through high lati- 
tude Galactic sight lines are also larger (124±8 K; 
Gillmon et al. 2006, or ranging from 81 K at 
log 7V(H2) = 20 to 219 K at log A^(H2) = 14; 
Wakker 2006). 

The population of J=3 rotational level is in 
turn enhanced compared to the Boltzmann dis- 
tribution, which indicates the presence of addi- 
tional excitation processes such as UV pumping 
(e.g. Noterdaeme et al. 2007a) and/or turbulent 
dissipation (as possibly indicated by the larger b- 
parameters for higher- J levels seen by Noterdaeme 
et al. 2007b). 

4.2. Density 

Absorption lines produced by neutral carbon 
are seen only in one component at Zabs = 3.28746 
(see Figure^]) and are associated with the strongest 
H2 component. 

We used the relative populations of the two first 
sub-levels of the C i ground state to derive the ex- 
citation temperature of the C i fine-structure level, 
according to the Boltzmann equation (see Eq. 1). 
We have adopted the energy difference between 
the C I excited (C i*: 2s^2p^ ^Pl) and true ground- 
state (2s22p2 3p0) levels, AE^g = 23.6 K. The 
population ratio 7V(J=1)/A^(J=0) of the Ci fine- 
structure level corresponds to an excitation tem- 
perature of Tox = 15.4±0.1 K. This is higher than 
the temperature expected in the case the excita- 
tion is dominated by the cosmic microwave back- 
ground radiation (Tcmbr = 11-7 K at z = 3.287) 
indicating that excitation by collisions is impor- 
tant. Using the results shown in fig. 12 of Srianand 
et al. (2005), we derive that the particle density, 
riH, is in the range 50-80 cm~'^. 

5. Discussion 

We have presented a detailed analysis of a 
QSO-DLA system with an extremely large col- 
umn density, log7V(Hi) = 22.0 ± 0.10, at Zabs = 
3.287 towards the quasar SDSS J 081634-^144612. 
The velocity structure of associated metal ab- 
sorption lines indicates the presence of 8 compo- 
nents grouped into three sub-systems centered at 
z = 3.28661, 3.28735 and 3.28814 respectively, 
spanning ~ 115 km s^^. Ci is detected in the 
z = 3.28735 sub-system whilst H2 is detected in 



both the z = 3.28661 and 3.28735 sub-systems. 
From the H2 excitation, we derive a kinetic tem- 
perature of Tk ~ 75 K and the observed column 
density ratio (C i)*/A^(C i) yields a particle den- 
sity in the range nn ~ 50 — 80 cm~'^. The deple- 
tion of metals onto dust grains measured in the 
strongest H2 component located at Zabs = 3.28735 
is similar to what is observed in the disc of the 
Galaxy. All this shows that this system, apart 
from having unusually large N(Hi), has proper- 
ties consistent with that of a typical H2-bearing 
DLA (see Ledoux et al. 2003; Noterdaeme et al. 
2008). 

While log A^(Hi) > 22 DLAs are very rarely 
seen in front of QSOs. several have already been 
detected in the optical afterglow spectrum and at 
the redshift of long-duration 7-ray bursts. Two 
of them (towards GRB 050401, GRB 080607 and 
GRB 060926) even have log A^(Hi) > 22.5 (Wat- 
son et al. 2006; Prochaska et al. 2009; Jakobsson 
et al. 2006). It is not surprising to observe a 
strong DLA at the redshift of GRBs since these 
objects are expected to be associated with star- 
forming regions where the gas is likely to be found 
in large quantities. Although QSO-DLAs are lo- 
cated close to regions where stars form as indicated 
by the presence of metals, the detection of C 11* ab- 
sorption and galaxy-like kinematics (Wolfe et al. 
2003), their exact nature is not completely eluci- 
dated. Some should be associated with the ISM 
of galaxies especially when molecules are detected 
(Noterdaeme et al. 2008), others are probably lo- 
cated in the outskirts of galactic haloes (M0ller 
et al. 2002; M0ller et al. 2004; Fox et al. 2007; 
Pontzen et al. 2008, Ranch et al. 2008; Rahmani 
et al. 2010; Fynbo et al. 2010, 2011). 

The difference between the populations of 
GRB-DLAs and QSO-DLAs is apparent because 
(i) the mean Hi column density is higher in 
GRBs than in DLAs reaching easily well beyond 
2Q21.5 Qj^-^ former case (see, e.g., Jakobsson 

et al. 2006; Fynbo et al. 2009) whilst QSO-DLA 
with log A^(H i) > 22 are very rare (Noterdaeme et 
al. 2009b), and (ii) the mean metallicity is larger 
for GRB-DLAs (Savaglio et al. 2006; Fynbo et 
al. 2006; Prochaska et al. 2007; Fynbo et al. 
2008), about 0.1 solar at z > 2 (to be compared 
to ~ 0.03 solar for intervening DLAs). Note that 
molecules (H2, CO) have been detected in a num- 
ber of QSO-DLAs (Noterdaeme et al. 2008, 2011, 
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Srianand et al. 2008), whereas they are rarely 
seen in GRB-DLAs (Ledoux et al. 2009). This 
may be a consequence of small number statistics 
however (see Fynbo et al. 2006 and Prochaska et 
al. 2009) and possibly of inadequate data as high 
spectral resolution in the blue is usually needed 
(see Ledoux et al. 2009). 

If these differences exist when comparing the 
overall populations, they may not be that appar- 
ent if we restrict ourselves to high H i column den- 
sity systems. This is why it is interesting to com- 
pare the properties of SDSS J 081634-1-144612 with 
those of GRB-DLAs. In Fig. [71 we plot the loga- 
rithm of the H I column density versus metallicity 
for QSO-DLAs where II2 is not detected (open 
circles), QSO-DLAs where II2 is detected (filled 
circles), the same for GRB-DLAs (squares) and 
SDSS J 0816344-144612 (filled diamond). It can 
be seen that there is a lack of systems with both a 
high metallicity and a high column density. This 
is well known for QSO-DLAs and could be a conse- 
quence of these DLAs being missed because of the 
high induced attenuation which makes the QSO 
drop out of the sample (e.g. Boisse et al. 1998). 
Altough GRB afterglows arc for a little while much 
brighter than quasars, the presence of a dust-bias 
could also affect their statistics (see e.g. Fynbo et 
al. 2009; Ledoux et al. 2009; and Greiner et al. 
2011). A possibility is that GRB-DLA mctallici- 
tics could be higher than measured as often only 
lower limits are derived from intermediate resolu- 
tion observations (see Prochaska 2006, Petitjean 
& Vergani, 2011). This can however probably 
not explain this lack of systems completely. SDSS 
J 081634-^144612 is located at the limit of the re- 
gion where systems arc missing and within the 
region where GRB-DLAs are located. In Fig. [H 
the logarithm of the H i column density is plotted 
for the same systems versus the depletion of iron 
onto dust grains. We have scaled the depletion 
of chromium in SDSS J 081634-M44612. It can 
be seen here again that the DLA towards SDSS 
J 081634-1-144612 is well within the region where 
GRB-DLAs are located. Therefore, our study in- 
dicates that extremely large Hi column density 
DLAs found toward QSOs have properties similar 
to those of GRB-DLAs. 

Following Schaye (2001), Krumholz et al (2009) 
proposed a radiation-insensitive model -hence ap- 
plying to both QSO and GRB-DLAs- where the 



lack of high-7V(Hi), high-metallicity is explained 
by the conversion from atomic to molecular gas. 
Here, the physical conditions in the _ff2-bcaring 
cloud indicate that we are still observing diffuse 
gas. This is not in tension with the above model 
since dense, cold and molecular gas is expected 
to have a small cross section and is not easily 
intercepted by the line-of-sight. Interestingly, a 
fully molecular cloud has been observed in the 
case of the logA^(Hi) = 22.7 DLA associated to 
GRB080607 (Prochaska et al. 2009), which could 
be observed very quickly (<lh) after the burst. 
Therefore, while the absorption properties of high- 
iV(Hi) GRB and QSO-DLAs appear to be similar, 
detecting fully molecular clouds will remain chal- 
lenging in the case of QSOs, because of the ran- 
dom distribution of the lines-of-sight and the high 
induced extinction. 

Very little is known about GRB host galaxies 
at high redshift (z > 2, Savaglio et al. 2009). 
Cosmological simulations show that GRB-DLAs 
are predominantly associated with haloes of mass 
10^° < Mvir/Mg < 10^2, an order of magnitude 
larger than the galaxies responsible for the bulk 
of QSO-DLAs (Pontzen et al. 2010; see however 
Barnes & Haehnelt 2010). But what is true for the 
overall QSO-DLA population does not hold for its 
high column density end. We have shown here 
that DLAs with the highest Hi column densities, 
seem to have absorption properties similar to that 
of GRB-DLAs. It would be most interesting to 
search for emission lines associated to these DLAs 
which have probably small inipact parameters (see 
e.g. Noterdaeme et al. 2012jl 

Note that the above conclusion is probably not 
very surprising as high H i column densities likely 
arise in gas located in the inner regions of star 
forming galaxies of moderate attenuation. This 
idea should be investigated with larger samples of 



•^The galaxy responsible for the recently found log7V(Hl) ~ 
22 DLA towards J 1135-0010 was detected at 6 ft! 0.1 are- 
sec from the background QSO (Noterdaeme et al. 2012). 
In this particular case, Kukarni et al. (2012) demonstrated 
that detecting the Lyman-a emission is even possible from 
UVES data. Here, the DLA is covered at the edge of two 
UVES echelle orders and possibly contaminated by scat- 
tered light in the red arm, preventing us from putting any 
meaningful limit on the Ly-a flux —within the slit. Follow- 
up observations specifically tuned to the search of emission 
lines is thus desirable (see e.g. Fynbo et al. 2010, 2011, 
Peroux et al. 2012). 
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both QSO and GRB-DLAs with high column den- 
sities. The BOSS survey (Eisenstein et ah 2011) 
will soon increase the number of such QSO-DLAs 
by an order of magnitude. 
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Table 1: Results of Voigt profile fitting analysis. 
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16.31+0.01 
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14.54+0.01 


13.53+0.01 
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radial velocity (km/s) 

Fig. 1.— Hi Lyman-a profile of the DLA at z = 3.286 toward QSO SDSS J 081634+144612. The overplotted 
solid line and accompanying dashed lines correspond to the best fit solution log Af(Hi) = 22.0±0.10. The 
origin of the velocity scale is taken at z = 3.287. 

Table 2: Molecular hydrogen column densities and excitation temperatures 
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1 With respect to Ci. 

^ The errors on the H2 column densities correspond to the best fits using the range of 6- values. 

3 Total log 7V(H2)= 17.60. 

4 Total log iV(H2)= 18.62. 
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Fig. 2. — Observed absorption profiles of Ci, Ci*, Znii, Fell, Cm, Mil, Sin together with fitted profiles. 
The dashed vertical lines indicate the positions of the velocity components. The origin of the velocity scale 
is taken at z = 3.287456, the redshift of the C i component. Regions affected by blends (sec Fell lines as 
well as the reddest clump of Cr iiA2066) arc not considered in the fitting process. Additional absorption seen 
at +150 kms^i the Ci* A1560.6 panel and in the red wing of CiA1656 and Ci* A1657.9 are due to Ci**. 
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Fig. 3. — Depletion pattern of elements relative to zinc measured in individual components in the DLA 
toward SDSS J 081634+144612. The points are slightly shifted along the x-axis according to their position 
in the system (for clarity purpose only). The filled circle indicates the component associated with the 
strongest H2 absorption. The dotted, dashed and continuous lines represent the typical relative abundances 
observed in, respectively, cold or warm gas in the Galactic disk and diffuse gas in the Galactic Halo, from 
Welty et al. (1999). 
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Fig. 4. — Left: The SDSS spectrum of J 081634+144612 (black) is shown together with the unreddcned 
SDSS composite (grey) and the same reddened using a SMC extinction law and E(B-V) = 0.05 (red). 
Right: Distribution of E(B-V) measured for a control sample of 163 QSOs within Az = ±0.02 from SDSS 
J 081634+144612. The arrow indicate the position of the latter. 
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Fig. 5. — Portions of the SDSS J 081634+144612 normalized spectrum. The best fit to the H2 absorption lines 
in several Lyman bands (0-0 to 4-0) is superimposed onto the spectrum, with the corresponding uncertainty 
represented by the shaded area. The fit parameters are given in Table [2] The blue and green labels indicate 
the branch and the rotational level for the bluest and reddest component of H2 , respectively. 
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Fig. 6. — H2 excitation diagram for the 2 components observed toward SDSS J 081634+144612. Triangles 
correspond to Zabs = 3.28667 and squares to Zabs = 3.28742. The column density Nj divided by the statistical 
weight, gj, is plotted for the J = up to J = 3 H2 rotational levels on a logarithmic scale against the excitation 
energy, Ej, in K. The J=0 to J=2 points have been fitted with straight lines using Eq.[T]giving the excitation 
temperatures T^^^ indicated in Table [2j 
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Fig. 7. — Logarithm of the total neutral hydrogen column density versus metallicity. Circles correspond to 
QSO-DLAs (Noterdaeme ct al. 2008), squares to GRB-DLAs (Prochaska et al. 2009), triangles to GRB- 
DLAs(Ledoux et al. 2009), inverted triangle to the QSO-DLA towards SDSS J 113520+001053 (Noterdame 
et al. 2012) and the diamond indicates our measurement for the QSO-DLA towards SDSS J 081634-1-144612. 
Filled symbols indicate systems in which H2 is detected, the filled square corresponds to the GRB080607 
(Prochaska et al. 2009). 
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Fig. 8. — Logarithm of the total neutral hydrogen cohimn density versus depletion. Circles correspond to 
QSO-DLAs (Noterdaeme ct al. 2008), squares to GRB-DLAs (Prochaska et al. 2009), triangles to GRB- 
DLAs (Ledoux et al. 2009), inverted triangle to the QSO-DLA towards SDSS J 113520+001053 (Noterdame 
et al. 2012) and the diamond indicates our measurement for the QSO-DLA towards SDSS J 081634+144612. 
Filled symbols indicate systems in which H2 is detected, the filled square corresponds to the GRB080607 
(Prochaska et al. 2009) . 
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